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Abstract 

We discuss the cosmology of string models with perturbative supersymmetry break- 
ing at a scale of O(TeV). Such models exhibit Kaluza-Klein like spectra and contain 
unstable massive gravitinos/gravitons. We find that considerations of primordial 
nucleosynthesis constrain the maximum temperature following inflation to be not 
much larger than the supersymmetry breaking scale. This imposes conflicting re- 
quirements on the scalar field driving inflation, making it rather difficult to construct 
a consistent cosmological history for such models. 
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There has recently been much interest in the possibility of realistic string theories 
with spontaneous perturbative super symmetry breaking |], |2[. These theories are very 
predictive in that they yield, in addition to the supersymmetric standard model, an en- 
tirely new phenomenon with a striking experimental signature |IJ-@. Specifically they 
contain a repeating spectrum of Kaluza-Klein (KK) modes all the way up to the Planck 
scale, whose spacing (e pa 1/2R, where R is the radius of compactification) is comparable 
to the supersymmetry breaking scale which is of O(TeV). These modes can be excited 
at forthcoming accelerators such as the LHC |fl, hence such models should be of im- 
mediate interest to experimentalists. However there are open questions concerning the 
cosmological viability of these models which need to be addressed first. In this paper we 
will investigate whether the Kaluza-Klein spectrum is consistent with cosmological con- 
straints on massive unstable relic particles ||. (Note that the more commonly discussed 
models with dynamical supersymmetry breaking in a 'hidden' sector are also constrained 
by similar cosmological considerations [7j §.) 

The only previous relevant discussion on the cosmology of Kaluza-Klein theories con- 
centrated on the prospect that they may include absolutely stable massive particles re- 
ferred to as 'pyrgons' 0. Such particles reside on the first rung of the ladder of Kaluza- 
Klein states, and are unable to decay because they carry a charge which is not exhibited 
by any of the massless particles. We shall not consider such models since, as we shall see, 
'pyrgons' do not exist in string theories with spontaneously broken supersymmetry. Al- 
though motivated primarily by string theory, our discussion will apply to all Kaluza-Klein 
theories in which the couplings are approximately independent of winding number. 

For such a theory, the thermal history of the universe is radically altered in the follow- 
ing way. We assume, as is usual, that there was an inflationary DeSitter phase, followed 
by reheating to a temperature Tr [lOfl. In conventional supergravity, reheating results 



in the production of gravitinos with number density proportional to T R ; the subsequent 
decays of the gravitinos can adversely affect primordial nucleosynthesis and requiring that 
they not do so results in an upper bound on the reheat temperature of ~ 10 5 TeV [[11] - 
[T3| . After reheating, the entropy, which we shall assume is subsequently conserved, is 
evenly spread out amongst the strongly (as opposed to gravitationally) interacting KK 
modes and the massless matter multiplets. At a temperature much higher than the KK 
level-spacing (T ^> e), nearly all the entropy is in the KK modes and almost none in the 
matter multiplets. Until the temperature drops below the first KK level, the evolution of 
the universe is therefore governed by the KK modes, whose contribution to the entropy is 
continually decreasing as the temperature drops. During this period there is production 
of massive gravitons and gravitinos which can only decay to the massless (twisted) parti- 
cles since their decays to untwisted KK modes is kinematically suppressed. Under these 
circumstances one might suspect that there is very severe bound on Tr and this indeed 
turns out to be the case. 

Let us first present the 'conventional' picture. Gravitinos are generated at high tem- 
peratures by two-body scatterings and the equation governing their number density is 

h 3/2 + 3H n 3/2 = (crv)n 2 Tad - , (1) 

T3/2 

where H is the Hubble expansion rate and (av) is the thermally averaged cross-section for 
gravitino production in the radiation bath of number density n rac j. The gravitino lifetime 



2 



r 3 /2 is given at rest by ]TJ 



,5 _ /™3/2\ 1 



r 3/2 ~ M^/m^ /2 ~ 10 5 sec , (2) 

where Mpi = G^ 2 ~ 1.22 x 10 16 TeV. Given the effective N = 1 supergravity couplings, 



27T — 

+\- ^d^Y^ + h.c. , (3) 



it can be shown that (av) ~ (Sn/Mp^ at temperatures T <C Mp\ [12]. The radiation 
density is given by 

n ra d = fl-(T) 5—, (4) 

where g(T) counts the relativistic degrees of freedom contributing to the total number 
density and is constant above temperatures of 0(TeV) for the minimal supersymmetric 
standard model (MSSM), with 

0(T a ) = 9 = 427/2 . (5) 

Assuming the canonical radiation-dominated evolution at this time, we can solve eq. ([TJ) 
to obtain 

v , n 3 / 2 g„(T) n I3d (T R )(av) 

Y 3/2 (T) = —!- = exp -£/r 3/2 ) , 6 

where time is related to temperature as 

t = 2.42 x 1(T 12 \g s (T)]- 112 sec , (7) 

and the factor g s (T)/g s (T R ) takes into account the decrease in the number of relativistic 
degrees of freedom, given constant total entropy 

tt 2 T 3 

si? 3 = q s (T) R 3 = constant . (8) 

y 30 V ' 

Note that we have taken g s {T) to also be the number of degrees of freedom determining 
the total energy density, as is appropriate at temperatures above a few MeV (when the 
neutrinos decouple). For the MSSM, one has 

g s {Tn) = § s = 915/4 , (9) 

at high temperatures when all particles are relativistic. 

Now let us consider the cosmological evolution when KK modes are present. Above 
the supersymmetry breaking scale the number of relativistic degrees of freedom is now 
no longer constant. The KK modes are labelled by quantum numbers of internal mo- 
menta/charges which are of the form 

n mR 

Vj^ — • (10) 
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where R represents some internal radius of compactiflcation. The winding modes (m 7^ 0) 
have masses of 0(Mp\) and need not be considered further, while the particles in the nth 
KK mode have masses m n ~ ne. Roughly speaking, whenever the temperature is raised 
by e, two new levels of (gauge interacting) KK excitations becomes relativistic, so that 
the number of degrees of freedom increases linearly with temperature. We can allow for 
this by writing 

T T 
g(T) = (g+ -g K ) , 9s(t) = (g s + -Qsk) ■ (11) 
e e 

The constants gK, g S K are determined by evaluating the number density and entropy 

density, respectively, of the plasma. For example consider the number density of KK 

modes in equilibrium: 

neq(T) = ? jtf I ^ exp(Vp 2 + n^/T)±i ' (12) 

where gi is the total number of interacting degrees of freedom in any KK level. Using 
various redefinitions, this becomes 

g{T A r°° , f°° , x 2 . , 

n eq (T) — / dx dy , (13) 

7r z e Jo Je/T exp y/X z + y z ± 1 

where we have approximated the sum at small e/T by an integral, and included a factor 
of two for positive and negative values of internal momentum (defined above as n). The 
integral becomes temperature independent for T ^> e, and we find that it deviates by less 
than 5% from the T 4 behaviour for high values of T/e (3> 2). The limiting values (when 
e/T — > 0) may be determined analytically to be 

9k = gi(XF + Xb) , (14) 

where 

Xf = 7vr 5 /480 C(3) = 3.71 , X b = vr 5 /60 C(3) = 4.24 . (15) 

In a similar fashion, the contribution of the KK modes to the total entropy may be used 
to determine g S K- Using s = (p +p)/T, we find that 

aT 4 f 00 , f 00 , x 2 (4x 2 /3 + 3y 2 ) , . 

s eq (T) « / dx / ^ 1 7=^= — I ^ (16) 

n z e Jo Je/T expv^ + y ±1 + y 



which gives limiting values of 



= £?i(XsF + Xsb) , (17) 



where 

Xf = 10125 C(5)/64vr 3 = 5.29 , xb = 675 C(5)/4tt 3 = 5.64 . (18) 

In the spontaneously broken string theories, each KK level comes in N = 4 multiplets, 
so that KK gauge bosons contribute 8 bosonic and 8 fermionic degrees of freedom in the 
vector and fermionic representations of SO (8) respectively. In the minimal case in which 
the KK excitations are in SU(3) <S> SU(3) C multiplets this gives 

g K = 128(xf + Xb) = 1018, g sK = 128( Xs f + Xsb) = 1400 . (19) 
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Of course there are additional contributions from higgs multiplets which are also expected 
to have KK excitations, so we shall consider these values to be a lower limit. To find the 
time-temperature relation, we assume that after inflation the metric is of the usual FRW 
form, with all the relativistic degrees of freedom in chemical equilibrium and therefore at 
the same temperature. Entropy conservation then gives 

SR 3 = ^K(T)7T 2 r 4 ^ 3 = constant; (2Q) 
6 oU 

and in particular 



4T 

H = —^- (21) 



3T 



Thus the Hubble parameter is 



T b/2 



Htr, = lM "-fjfc- 1221 



Differentiating with respect to time and substituting eq.(^), we find 



t(T) = — — . (23) 
15H(T) 

For the minimal value of g S K above, this becomes 

/ , s 1/2 , J> . -5/2 

t = 6.9xl0-"sec(— ) (— ) , (24) 

at temperatures T ^> e. 

In order to ascertain the abundance of massive gravitons/gravitinos, we need to iden- 
tify the processes which can contribute to their manufacture and decay. Vertices between 
KK modes (which come from untwisted sectors of the string theory) must satisfy the con- 
dition that their internal moment a/ charges are conserved, and are simply proportional 
to the string coupling constant. In addition vertices can exist between untwisted modes 
and twisted (massless) matter multiplets, with couplings g n oc g<5~ m «/ M pi m g, where 5 
is some constant depending on the type of compactification [nj. For masses much less 



than Mpi these are clearly unsuppressed, so that we can write down effective terms for 
the coupling of the KK modes to each other, and to the massless multiplets. The creation 
of KK gravitons and gravitinos goes via effective four particle interactions. For example 
the cross section for 

A a k + A\ - X c m + V n (25) 

goes as 

o ~ ^f abc f abc ^(n + m -k-l), (26) 

647T Mpj 

where the integers k,l,m,n label KK modes, and we have omitted numerical factors of 
(9(1) coming from the trace over solutions to the Rarita-Schwinger equation and phase- 
space integrations. We therefore write the total cross section for n-gravitino production 
from k plus I interactions as 

a k i mn = a 5(n + m - k - I) . (27) 
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where a is a factor of 0(8ir/M^{), which is dependent on the details of the model, but in- 
dependent of the KK-modes. In addition there is a contribution coming from the massless 
sector which we shall neglect since there is only one such sector. 

The massive n-gravitinos/gravitons may decay into either two massless twisted states, 
or to untwisted I plus n — /-states. In the first case the decay rate is found to be, 



twisted 



(28) 



The decay to untwisted states is kinematically suppressed however. Consider a positive-n 
state decaying to an /-state plus an n — /-state. If / is negative then the masses of the 
products is \n\e + 2|/|e which is larger than the mass, \n\e, of the decaying particle. If / 
is also positive, then sum of the product masses is equal to that of the decaying particle. 
The decay rate is proportional to the momentum of the decay products in the centre of 
mass frame: 



r„„_ « H = Ji + " 2 - ' V ') 2 ' 2 - 2 E±<^m iM = o . (29) 

V n 4 n z 2 

As discussed in ref. ||, there may be a mass splitting of 0(e) in any KK level, so that the 
decay rate to untwisted modes will generally be suppressed by a factor l/\n\. When the 
sum over / is taken this decay mode may be of the same order as the decay to twisted 
states. We shall therefore take the lifetime of n-gravitons/gravitinos to be 

r n ~9.8x 10 4 sec(^) ' |n|" 3 , (30) 
corresponding to the lifetime for decay of standard gravitinos into photons and photinos 



|14| . Inclusion of all the strongly interacting, twisted, final states would speed up the 
decays by a small factor. For example, if the twisted products consisted of all the matter 
and Higgs particles in the MSSM, then the above would be reduced by a factor 12/49. 

With these estimates we are ready to tackle the evolution of the n-graviton/gravitino 



number density n n . This is governed by the equation ]lf 

, o 77 , n n r r r r d3 ^ d3 qi d3( ^ d3 q^ t , \ 

„„ + 3Hnn + - = e e y y y y m , ^ ^ ^ w * <* +«-*.-*) 

fkfl(l — fm) \Mkl^ mn \ 2 (31) 

where the fi are the occupation numbers and we have assumed that the massive particles 
are non-relativistic when they decay. This equation describes the creation of a gravita- 
tionally interacting n-state plus a gauge interacting m-state, from gauge interacting k 
plus /-states. We have taken the occupation number of the n-state to be negligible (since 
it is never in equilibrium), so that the reverse process does not occur. We may reasonably 
adopt equilibrium distributions for the three remaining (k,l,m) states, omit the Pauli 
blocking factor for the m-state, and rewrite eq.([H]) as 



h n + 3Hn n + — = a kimn n eq (k,T) n eq (i,T), (32) 

' " mkl 
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where, since we are concerned with the longest lived states, we have taken v « c. Although 
approximate, this expression has the correct temperature dependence and we shall bury 
our approximations in the parameter a which is still of 0(87r/Mpj). Defining 

n n = n n exp(-t/r n ), (33) 



and using eq. (|21|) , we find 




= ^<?kimn n e q(k,T) n eq (i,T) exp(-£/r n ) . (34) 

mkl 

The number density of the individual k and I states follow some equilibrium curve which 
obeys, e.g. 

+00 

™rad ~ n eq (T) = n eq(k,T) ■ (35) 

k=— 00 

Since the epoch of nucleosynthesis is much later than the time at which these particles 
are created, we take the exponential factor to be unity for the purposes of calculating 
the initial abundance. Performing the summations in eq.(^) and integrating from Tr to 
T = e gives 

lim = |^^W«p(-«/r.), (36) 
9 g s g H(T K ) 

where, again, there is a factor to account for the change in photon number density between 
e and the final temperature T < e. Note that the final density of n-gravitons/gravitinos is 
proportional to T^ 2 . For typical parameter values (a = 8ir/Mp h = 1018, g S K = 1400), 
this is 

/ , X 1/2 / rp s 3/2 

Y ni T) ~ 3 x 10- 16 (— ) exp(-t/r B ) (37) 

at T < e; the relic abundance during nucleosynthesis would be smaller by a factor g s (T)/g s (T ) 
where g s (To) = 43/11 is the effective number of entropic degrees of freedom at T <ti m e 
(taking into account the three decoupled relativistic neutrino species). This is close to 
the value of 3.36 for the effective number of degrees of freedom contributing to the total 
energy density, so for convenience we have ignored the small difference. 

The energy density in the relic KK modes decreases as T 17//3 for T > e, i.e. faster than 
the energy density in 'radiation' (including KK excitations) which goes as T 5 . (At T < e 
the former decreases as T 16//3 , while the latter does so as T 4 .) Therefore the universe will 
become 'radiation' dominated when 

7T 2 T 4 

\n\ e Y n (T) n rad (T) < g s {T) — — . (38) 

where g s {T) is taken from eq.(|TI|). Anticipating the bounds which we will find on Tr, the 
temperature at equality will be higher than e, therefore the appropriate time-temperature 
relationship is eq . fl24]) . Since at late times, the number density of KK modes is dominated 
by the lighter (slowly decaying) particles which have a small exponential suppression factor 
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in Y n , we may approximate the sum on the left above by an integral, 

t \ - 2 / 3 



— poo 

\n\ exp(— tv? /ti) »s / d 2 x x exp(— x 3 ) 

Jo 



n 



2 rtt)^r' 3 . (39) 



3 \3J Vri 

This expression is valid for t^Ti, at later times, the KK modes have all decayed, i.e. 
their number is exponentially suppressed. Therefore 'radiation' domination will occur at 

/ , \ 5/2 / T \ -1 

r>l.lxlO«TeV — ^ (40) 



corresponding to a time 



, e \i/2 / r D \ 5 / 2 
t <5.3 x 1(T 29 sec — - -^t • (41) 



TeW VTeV 

It may appear surprising that the higher the reheat temperature, the less likely the grav- 
itinos are to matter dominate at a late time. This can be explained as follows; there is 
only a limited amount of entropy, and when Tr is high, more of it is initially distributed 
in heavier modes with higher KK number. Since these modes decay more rapidly, they 
are able to matter dominate only at very early times. 

We can now examine the effect of the decaying particles on the abundances of the light 
elements. The effect of hadronic decays occuring at the beginning of the nucleosynthesis 
era on neutron-proton interconversions has been studied in detail in ref. |17]]- in the time- 
interval 1 <, t <, 10 2 sec, the requirement that the 4 He mass fraction not be increased above 
25% translates into the requirement (see figure 3 in ref. 0): 

X>| F„(T)n rad (T) <2xKT n , (42) 

which corresponds to the bound 

T R <21TeV^— J • (43) 

For t > 10 4 sec, the photodissociation of 2 H due to the radiation cascades triggered by the 
decaying particles impose the constraint 

E W Y nV) ^rad(T) < 1(T 13 , (44) 

n TeV 

or 

T R <2.7TeV(— J . (45) 

One should however bear in mind that the approximation we used to calculate the KK 
number density breaks down when Tr becomes comparable to e since the number density 
is then suppressed by a Boltzmann factor. We can however justifiably assert that the 
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maximum temperature which the universe reached cannot have significantly exceeded the 
supersymmetry breaking scale. 

A precise formulation of the early history of the universe is lacking for these models, 
however it is clearly essential that there be an inflationary phase ||10|| , followed, as we 



have shown, by reheating to a temperature no greater than the supersymmetry breaking 
scale. We now argue that this imposes a conflicting set of requirements on the scalar 
field $ which is presumed to drive inflation. In order to account for the amplitude of 
the scale- invariant density fluctuations probed by COBE \TE\, m$ is required to be of 
O(10 s ) TeV fpL9| . For the reheat temperature to be low, the inflaton must then be very 



weakly coupled to matter fields. However even assuming only gravitational couplings, i.e. 
a decay rate T$ ~ m|/Mp[, the reheat temperature cannot be reduced below 



R 



( r t M Pl) v^ 103 TeV ^ { _^_ } , (46) 



where Tr has been obtained by equating T$ to the Hubble rate H. In fact, the only 
mass scale in these models, apart from the Planck scale, is the supersymmetry breaking 
scale. One could then envisage a scenario with two epochs of inflation ||; the first stage to 
create the correct level of density fluctuations, and the second, with m$ ~ e, to remove the 
KK states. However assuming gravitational couplings as above, the reheat temperature 
is then of 0(keV), i.e. too low for even primordial nucleosynthesis to occur. On the 
other hand if the inflaton is coupled directly, the reheat temperature would be too high, 
unless the gauge coupling is unnaturally small, of O(10~ 4 ) ||. A way out may be to 
introduce an intermediate scale of unification, with a low-mass 'flaton' singlet remaining 
after symmetry breaking ||20|| . However the reheat temperature is then of 0(MeV), i.e. 



barely high enough for nucleosynthesis, and further, the Affleck-Dine mechanism must be 
invoked in order for the baryons to be synthesized just beforehand. 

In conclusion, it appears difficult to construct a consistent cosmological history for 
models with spontaneous supersymmetry breaking. In view of their many other advanta- 
geous features Q, possible resolutions to this problem should be pursued. 
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